Thiols in hydrothermal solution: Standard partial molal properties and 
their role in the organic geochemistry of hydrothermal environments 

Mitchell D. Schulte*'* and Karyn L. Rogers 5 

*NASA Ames Research Center, Exobiology Branch, Mail Stop 239-4, 

Moffett Field, CA 94035 

^Washington University, Department of Earth & Planetary Sciences, Campus Box 1 169, 
One Brookings Drive, St. Louis, MO 63130 


‘Author to whom correspondence should be addressed (mschulte@mail.arc.nasa.gov) 


Abstract - Modem seafloor hydrothermal systems are locations where great varieties of 
geochemistry occur due to the enormous disequilibrium between vent fluids and seawater. The 
disequilibrium geochemistry has been hypothesized to include reactions to synthesize organic 
compounds (Schulte and Shock, 1995; Shock and Schulte, 1998). Despite the incomplete 
understanding of the carbon budget in hydrothermal systems, the organic geochemistry of these 
sites has received little attention. Experimental simulations of these environments, however, 
indicate that organic compounds may have difficulty forming in a purely aqueous environment. On 
the other hand, thiols, thioesters and disulfides have been implicated as reaction intermediates 
between CO or C0 2 in experiments of carbon reduction in hydrothermal environments, as well as 
in a variety of biological processes and other abiotic reactions (Wachtershauser, 1990; Heinen and 
Lauwers, 1996; Huber and Wachtershauser, 1997; Russell et al., 1998). The reduction of C0 2 to 
thiols, for example, is observed using the FeS-H 2 S/FeS 2 couple to provide the reducing power 
(see Schoonen et al., 1999). We have used recent advances in theoretical geochemistry to estimate 
the standard partial molal thermodynamic properties and parameters for the revised Helgeson- 
Kirkham-Flowers equation of state for aqueous straight-chain alkyl thiols. With these data and 
parameters we have evaluated the role that organic sulfur compounds may play as reaction 
intermediates during organic compound synthesis. We conclude that organic sulfur compounds 
may hold the key to the organic chemistry leading to the origin of life in hydrothermal settings. 
These results may also explain the presence of sulfur in a number of biomolecules present in 
ancient thermophilic microorganisms. 


INTRODUCTION 

Sulfur may play and may have played a key role in abiotic organic synthesis and 
biochemistry in hydrothermal environments (de Duve, 1991; Russell et al., 1998), and recent 
experimental work to model these environments appears to support this notion (Heinen and 
Lauwers, 1996; Huber and Wachtershauser, 1997; Cody et al., 2000). Furthermore, the presence 
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of sulfur in a number of key biological molecules (such as the amino acids cysteine, methionine 
and cystine, and coenzymes M and acetyl CoA) and metabolic processes (sulfur oxidation and 
reduction in a number of chemolithoautotrophs) also lends credence to this idea. The large number 
of oxidation states possible for sulfur makes for diverse chemistry, especially in combination with 
carbon. Thiols, the simplest sulfur-bearing organic compounds, are also common in biochemical 
processes and are, for example, responsible for the distinctive odors of garlic ( 2 -propene-l -thiol), 

onions (propanethiol) and skunks ( 3 -methyl- 1-butanethiol). 

Understanding the biochemical processes of the thermophilic microorganisms that inhabit 

hydrothermal environments and that can live to temperatures at least as high as 1 13 ( 

al„ 1997) requires knowledge of the properties of the chemical reactions involved at elevated 
temperatures and pressures. In order to assess the role of aqueous organic sulfur compounds in 
hydrothermal organic geochemistry, we have been investigating their thermodynamic properties. 
We have culled the literature to obtain die experimentally deteimined properties of organic sulfur 
compounds. We were able to calculate a number of essential properties, such as free energies of 
formation, from solubility data available in the literature together with standard properties of 
organic sulfur gases. However, a number of the properties of aqueous organic sulfur compounds 
have not been experimentally determined. Furthermore, most of thermodynamic data that are 
available are only for 25°C and 1 bar. In order to determine reaction properties at die temperatures 
and pressures of the hydrothermal systems in which thermophilic organisms live, we use equations 
of state developed by Helgeson and co-workers (Helgeson et al., 1981), a summary of th 
equations is given in the Appendix. A key piece of intonation needed to extrapolate 
thermodynamic properties to elevated temperatures is the partial molal heat capacity, for which 
experimental data are unavailable for nearly all aqueous organic sulfur compounds, even at 
standard state conditions (25"C and 1 bar). We have used correlation methods to esdmate the 
partial molal heat capacities and volumes of many organic solutes. These estimates allow us to 
asses the role of organic sulfur compounds during the reduction of carbon in hydrothermal 

settings. 
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ESTIMATION METHODS 


Methods have been developed by Helgeson and co-workers (Schulte et al., 2001; Shock 
and Helgeson, 1988, 1990; Shock et al., 1989, 1992, 1997 ; Sveijensky et al., 1997) to estimate 
the standard partial molal thermodynamic properties of aqueous species. These methods have been 
extended to include aqueous organic compounds (Dale et al., 1997; Prapaipong et al., 1999, 
Schulte and Shock, 1993, 2001; Shock and Helgeson, 1990; Shock and Koretsky, 1993, 1995). 
Recent high temperature experimental measurements of the partial molal volumes and heat 
capacities of a selected group of compounds (Criss and Wood, 1996; Schulte et al., 1999) has led 
to improvements in these methods (Schulte, 1997; Schulte and Shock, 2001). Described below 
are the methods by which the standard partial molal thermodynamic properties of the aqueous alkyl 
thiols were calculated or estimated. 

AGf, AH° f and S 

Values of the standard partial molal Gibbs free energy of formation (AG^) for aqueous 
methanethiol (CH 3 SH) to butanethiol (CH 3 (CH 2 ) 3 SH) were calculated using the standard Gibbs 
free energy of formation of the gaseous species (Alberty et al., 1987) along with the Gibbs free 
energy of reaction for the hydration of the gas into the aqueous phase (Abraham et al., 1990). 
Similarly, the standard partial molal Gibbs free energy of formation for aqueous pentanethiol was 
calculated using data from Alberty et al. (1987) and Amoore and Buttery (1978). The standard 
partial molal enthalpies of formation (AH}) for aqueous methanethiol to butanethiol were 
calculated in a similar manner, using the standard enthalpies of formation of the gaseous species 
(Alberty et al., 1987) and enthalpies of hydration reactions (Abraham et al., 1990). Figures 1 and 
2 show the resulting values of the standard Gibbs free energy of formation and standard enthalpy 
of formation, respectively, for the aqueous thiols as a function of the number of carbons in the 
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aJkyl chain (i.e. 1 for methanethiol (CH s SH), 2 for ethanethiol (C^SH), etc.). As is evident in 
each figure, the standard partial molal thermodynamic properties of the alkyl thiols, like for all 
organic compound families, display systematic behavior as alkyl chain length increases. Thts 
behavior is due to the addition of the -CH 2 - group in increasing the chain length. The addition of 
this group increases the standard partial molal Gibbs free energy of formation of aqueous organic 
species by 2050 cal mol 4 and decreases the standard partial molal enthalpy by 5670 cal mol 
(Shock and Helgeson, 1990). The straight lines shown on each plot represent best fit lines to the 
calculated data (shown as symbols in each plot) using 2050 cal mol 4 as die slope for the standard 
partial molal Gibbs free energy of formation and -5670 cal mol 4 for die standard partial molal 
enthalpy of formation. The best line fits to die calculated date using these slopes results tn 
intercepts of the lines of 4370 cal mol 4 for AG} and -6180 cal mol 4 for AH,, allowing 
calculate these properties for suaight-chtened aqueous Alois for which experimental date are not 

available. The equations used to make these estimates are 


AG°= -4370 + 2050n 


( 1 ) 


and 


AH} = -6180 + 5670n , 


( 2 ) 


where n represents the number of carbons in the alkyl chain. 

values of the standard partial molal entropy (S’) of aqueous methanethiol to butanethiol 

were calculated using the values of AG] and AH}, along with values of S“ of the elements from 
Cox et al. (1989). Figure 3 shows the resulting values of the standard partial molal entropy as a 
function of the number of carbons in the alkyl chain, along with a best fit line to these values using 
the slope of 6.7 cal mol 4 K" taken from Shock and Helgeson (1990) for aqueous organic species. 

The line corresponds to the equation 
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S°= -33.1 + 6.7h, 


( 3 ) 


allowing us to calculate S° for aqueous thiols for which calculation of the entropy from the free 
energy and enthalpy, as described above, cannot be made. The values of S for methanethiol to 
butanethiol, obtained as described above, are listed in Table 1. 


C P and V 

Because there are no experimental data for the standard partial molal volumes and heat 
capacities of aqueous thiols, even at the standard conditions of 25°C and 1 bar, we have estimated 
these properties for the aqueous thiols. In order to estimate the heat capacities, we compared the 
standard partial molal heat capacity of the aqueous amino acid serine to that of the aqueous ammo 
acid cysteine. These two amino acids are identical in structure, except for the functional group at 
the terminal end of the alkyl chain). In serine, the terminal functional group is a hydroxyl (-OH). 
In cysteine, the terminal functional group is a thiol (-SH). The experimentally determined value of 
the standard partial molal heat capacity of serine is 28.1 cal mol 1 K' 1 (Jolicoeur and Boileau, 
1978), while the standard partial molal heat capacity of cysteine is 44.9 cal mol 1 K 1 (Jolicoeur et 
al., 1986). Therefore, the difference in the standard partial molal heat capacity of these two 
aqueous molecules is 16.8 cal mol 1 K 1 . This value is, then, the difference in standard partial 
molal heat capacity between a molecule containing a thiol group relative to that containing a 
hydroxyl group. Because the amino acids that were compared each contain three carbons in their 
alkyl chains, we have chosen to add the difference in standard partial molal heat capacity, 16.8 cal 


i The structure of serine is HO-CH 2 -(CH)NH 2 -COOH, while that of cysteine is HS-CH 2 

(CH)NH,-COOH. Note that the difference in the structures is the left-hand terminal functional 
group (-OH in serine and -SH in cysteine). 
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mol -t K-. ,o that of propanol to estimate the hea, capacity of ptopanethiol. The resulting value of 

the standard partial molal hea, capacity of aqueous propanethiol is 101.1 cal mol ' K". Values of 

C- for Ore other aqueous thiols at 25°C and 1 bar were then estimated by addmg or subtracting 
multiples of 21.2 cal mol ' K' from this value, following the arguments of Shock and Helgeson 
(1990) on die systematic behavior of partial molal hea, cavities of aqueous organ, c compoun 

families. 

Values of the standard partial molal volume for die aqueous alkyl thiols were eshmated ,n a 
manner similar to tha, used for die heat capacities. The difference in standard partial molal volumes 
between serine (V = 60,7 cm' mo, '; Cabani e, al„ 1981, and cysteine (V - 73.44 cm mob : 
Millero et al„ ,978) is ,2.87 cm' mo, '. This value was added to the partial mola, volume o^ 
propanol to esdmate the partial molal volume of propanethiol, resulting in a value of 83.57 cm 
mol- 1 The partial molal volumes of aqueous organic compounds also vary ^ ^ 

systematically wtih carbon number (Shock and Helgeson, 1990), and multiples of 15.8 cm' mol 
were added or subtract*, from this value as a taction of carbon number to estimate the standard 

partial molal volumes of other aqueous thiols at 25°C and 1 bar. 

To gain an appreciation of how accurate the estimated values may be, we have plotted the 

c; and V of the aqueous thiols against the number of carbons in the alky, chain in Figure 4. We 
have added to the plot the experimentally determined values for aqueous H 2 S (winch plots at n 
0). The values of <T P and V we would estimrue by the methods described above for n = 0 are 
37 5 cal K ' mof and 37.17 cm' mol ', respectively, while the experimental values for H,S are 
42 7 cal K ' mol ' and 34.92 cm' mol respectively. Our estimared values are consistent wtth the 
experimental values for hydrogen sulfide, giving us confidence in our estimation methods. 


Helgeson-Kirkham-Flowers equation of state parameters 

~Z , • Hfl rH pu _CH while that of propanethiol is HS-CH 2 - 

of propanol is HO-LH 2 -Ln 2 ^n 3 , , . i as t tv e terminal 

CHj-CHj. The difference in *^_™°“ta»unds f^onal group in 


The structure 


functional group in propanol, and the presence 
propanethiol. 
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The standard state thermodynamic properties estimated for aqueous alky lthiols at 25°C and 
1 bar listed in Table 2 were used to obtain parameters for the revised Helgeson-Kirkham-Flowers 
(HKF) equations of state (Helgeson et al., 1981; Tanger and Helgeson, 1988, Shock et al., 1992) 
using correlation algorithms for aqueous nonelectrolytes (including neutral organic species) 
presented by Shock and Helgeson (1990), Schulte et al. (2001) and Schulte and Shock (2001). 
The resulting parameters are also listed in Table 2. These data and parameters were used with the 
revised HKF equations of state (see appendix) to calculate standard state properties at elevated 
temperatures and pressures. Such calculations are facilitated by using the SUPCRT92 computer 
program (Johnson et al., 1992), which is available free of charge from Prof. Harold Helgeson at 

the University of California, Berkeley. 

FORMATION OF AQUEOUS THIOLS 
IN HYDROTHERMAL ENVIRONMENTS 

Surface promoted catalysis on sulfide minerals may lead to the formation of thiols 
hydrothermal environments (see e.g. Huber and Wachtershauser, 1997, Russell et al., 1998, y 
et al., 2000). The following reactions describe the formation of aqueous methanethiol from CO 
(Equation (4)), CO, (Equation (5)), H 2 S and H 2 , the most likely starting materials for organic 
synthesis in hydrothermal environments, especially on the early Earth. 

CO(aq) + H 2 S(aq) + 2 H 2 (aq) = CH 3 SH(aq) + H 2 0 

methanethiol 

C0 2 (aq) + H 2 S(aq) + 3 H 2 (aq) = CH 3 SH(aq) + 2 H 2 0 . ( 5 ) 

methanethiol 
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The carbon in CO or COj is reduced by the hydrogen and hydrogen sulfide that is p 
reactions among sulfur-bearing minerals (e.g. pyrite and pyrrhotite,. Urnger chain thiols can be 
formed similarly, according to the general reaction shown on Figure 5. Ftgure 5 depicts e 
logarithms of the equilibrium constants of the reactions for thiol formation from CO 
methanethiol (C.) to decanethio! (C, 0 ) as a function of temperature a. saturadon pressure. Figure 
shows these for thiol formation from CO, These plots indicate that die formauon of thiols rom 
simple starting materials is energedcally favored a, temperatures commonly encountered m 
hydrothermal environments. In addidon, i, is evident from comparing die curves in Ftgures an 

6 that reacdons that form thiols from CO are more diermodynamically favorable (log K is more 

■ i pci anf i jc favorable to higher 
positive) than the corresponding reactions starting 2 ’ 

f rn 900°C for formation from C0 2 ). This result is 
temperatures (~300°C for formation from CO vs. -200 C tor to 

consistent with .guments put forward by Schoonen e. al. (1999, in which the role of elec™, 

source of reduced carbon in hydrothenmal envirorunents. Ftgure 5 also indicates that format, on o 
longer chained compounds is more dienuodynamMly favomble dian diose of die shorter chaine 
compounds to a temperature of about 300“C. Similar conclusions were reached by Shoe an 
Schulte (1998) in their study of die potendal for organic synthesis dunng flutd nuxtng ,n 
hydrothermal systems. Based on these results, we have emulated dte potential for product, on o 
the simples, thiol, methanethiol. in modem and early Earth hydmthetmal systems, taktng tnto 
account die geochemical constraints imposed by modem and early Earth hydrothetma, 

environments. 


Modem Systems 


Geochemical constraints consist of naturally buffered aedvides and fttgacties of die 
chemical species that partietpate in these reactions. Measurements of die aqueous hydrogen 
concentrations of vent fluids are consistent with die pyrite-pyrrhotite-magnedte (PPM) turner 


Q 


assemblage (Shock et al , 1995). Therefore, model calculations for modem hyrothermal systems 
assume activities of H, and H,S in equilibrium with the PPM mineral assemblage (see Table 3). 
For the formation of methanethiol from CO, the activity of CO in modem hydrothermal systems is 
assumed to be a relatively consetvative value of 10 ‘ (log ^ = -6). For calculations involving 
CO,, the activity of CO, is set to 13.8 mM, taken from a compilation of measured CO, 
concentrations for modem seafloor hydrothermal vent systems (Von Damm, 1990). For the 
following model calculations, the activity of H,0 is assumed to be unity (Helgeson et al., 1993) 
and we have assumed a total pressure of 250 bars, consistent with a 2.5 km deep seafloor 

hydrothermal system. 

The final geochemical constraint we place on our model calculations is the equilibrium 
constant for the reaction. We use the relation 


log K = 


ag; 

-2.303RT 


( 6 ) 


obtain values of the logarithm of tire equilibrium constant (K) a. appropriate temperatures and 
pressures, where AO; is the standard Gibbs free energy of the reaction a. elevated temperature and 
pressure, R is the gas constant and T is die iemperature in Kelvin. For each reaction the law of 
mass action is used to determine the equilibrium activity of methanethiol. For example, at 

equilibrium 


( a CH 3 SH ) ( a H 2 o) 

(7) 

( a Co)( a H,s)( a H, ) 



where K, stands for the equilibrium constant of reaction (4) above and a, refers to the activity of 
each chemical species (i) involved in tire reaction, raised to the stoichiometric reaction coefficient 
(which is positive for products and negative for reactants). Taking the logarithm and rearranging 
the law of mass action of reaction (4) gives 
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log a CH SH = logK, + log a c0 + loga„ !S + 21oga„, , 


( 8 ) 


which allows us to assess the activity of aqueous methanethiol given the geochetrucal consents 
discussed above. Similarly, we have also calcuiated the equilibrium ntethanethtol achvthes or 
fonnation b, CO, as illustrated by Equation (5). The corresponding law of tnass acnon reduces to 


log a ( 


ch 3 sh 


= log K 5 + log a c02 + log a HjS + 3 log a Hj 


(9) 


The logarithm of dte aqueous tnethanethio, acdvity can he calculated using Bquahon (9, and dte 
geochemical constraints discussed above. 

The potential for ntethanethiol formation from CO and CO. can be calculated as a frncnon 
of temperature using Equations (8) and (9) and values of the geochemical constraints adjusted for 
temperature. The tesulting values of dte logaridtn, of dte equiiibrium methanethtol acnvtty om 
teacdons (4) md (5) are shown as a frtnctton of temperature in Figure 7. Examination of Bgure 

reveals that if methanethiol forms from CO, its acnvtty can be as high as 
temperatures), while fonnation of methanethiol from CO. gives an activity as high as 1 . o e 

that the calculated equilibrium acdvides of nterhanedtio, am much greater from formtng by 
reduction of CO than CO. by 4 ordetx of magnitude. The highest calculated acdvides are at dte 
iowest temperatures, md fdnetic bakers may prevent the fonnation reactions from prtteee tng 
under dtese conditions. In addition, dte calculated acdvides assume that dte CO or . urn y 
converted into methanethiol. and does no, constder the formation of other carbon compounds or 
consumption of CO and CO. in the system. Therefore it is unlihely ft* dtese high concennattons 
would be reached. On the other hand, the catalytic activity of the sulfide miner 
these systems may help overcome dtese barriers and stgniftcant amounrs of thiols may m ee orm 

and accumulate in these environments. 
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Archean Systems 


Hydrothermal systems were likely to be present on the early Earth, and may have 
more prevalent than in the modem day. Furthermore, the rocks controlling the composition of the 
vent fluids may have been more reducing than modem mid-ocean ridge basalts (MORBs). This 
would result in more hydrogen production during alteration of the crust, and the potential for even 
greater organic synthesis (Schulte and Shock. 1995). Therefore, we have performed calculations 
similar to those above using a set of geochemical constraints appropriate for Archean hydrothermal 
systems. We assume that the hydrogen activity in ancient hydmthermal systems was buffered by 
the FMQ (fayalite-quartz-magnetite) mineral assemblage (see Table 3). Since the FMQ min 
buffer does not itself contain any sulfur-bearing minerals, we must couple the FMQ buffer 
pyrrhotite-magnetite (PM) mineral assemblage to constrain the aqueous H,S activity. We use the 
activity of aqueous H, required by FMQ equilibrium as a function of temperature from 0-350 C to 
determine the aqueous H,S activity in equilibrium with the pyrrhotite-magnetite reaction (see Table 
3) as a function of temperature. The values of CO and CO, activity are assumed to be consistent 
with models of increased amounts of these gases in die Archean atmosphere. It has been 
postulated that the partial pressures of CO and CO, were as high as 1 and 10 bars, respectively, on 
the early Earth (Kasting. 1993). We have adopted these values for simplification for our 
calculations of the methanethiol activity in Archean hydrothermal systems. Using the geochemical 
constraints of Archean hydrothetmal systems described above, along with Equations (8) and (9), 
we can calculate the equilibrium methanethiol activities as a function of temperature. The results 
are shown in Figure 8, Because of the greater reduction potential in the early Earth systems, the 
calculated equilibrium activity of methanethiol is much higher than in the case of modem systems. 
Formation of methanethiol from CO gives activities as high as 10 s , while formation from CO, 
gives activities as high as 1C 1 . Again, low temperatures may prevent these reactions from reach, ng 
equilibrium; however, methanethiol activities at higher temperature are still significant. The 
potential is much greater to form organic compounds from CO than from CO.. 
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THIOLS AS REACTION INTERMEDIATES IN MODERN 
AND EARLY EARTH HYDROTHERMAL SYSTEMS 


Recent experimental and theoretical evidence indicates that thiols may be necessary 
intermediates for the formation of a q neous organic molecules. snch as carboxylic acids, ammo 
acids, hydrocarbons and .cohols, in hydrothermal envimnments (Heinen and Lauwers, * * 
Huber and Wachtershauser. 1997; Russell et al„ 1998; Schoonen et al„ 1999; Cody et al„ 

The formation of these molecules may involve the attachment of precursors (such as CO, CO, or 
formaldehyde) to sulfide mineral surfaces, followed by the reduction of the attached molecule to 
thiols, then formation of carboxylic acids, alcohols, or amino acids horn the intermediate thiols 

In order to assess the role of aqueous thiols as intermedia* species, we have calculared the 
equilibrium constants of some reactions involving simple aqueons thiols, alcohols and carboxylic 
acids in hydrothermal environments. The reactions forming methano, Dorn methanethio, an 
ethanol and acetic acid from ethanethio, are shown in Table 4. The equilibrium constants for these 
reactions are p.otred as a function of temperarure a, 250 bare in Figure 9. Note that in each 
formation of the products is favored at all temperatures relevant to hydrothermal systems. s 
indicares that once formed, thiols would react to produce other organic species. If we now 
into account die s„e geochemical consents discussed above for duel formation, we can 
calculate die amount of organic species diat would form in hydrothermal systems from die 
previously calculated activities of thiols. 

As shown previously, the equilibrium thiol acdvity depends on the geochemical consents 
imposed by die system and the carbon source for thin, fonnadon (e.g. CO or CO,,. In order to 
assess die production of methano, from methanethio, in modem and early Hard, hydrotherm 
systems, we used the equilibrium methanethio, acdvides calculared previously and shown m 
Figures 7 and 8, along with geochemiC consents imposed by die environment. In addition, we 
calculated equilibrium activities of ethanethiol for formation from CO and C0 2 in 





the production of ethanol and acetic acid from ethanethiol. The results for both modern and early 
Earth hydrothermal systems, with the formation of the intermediate thiol from both CO and CO,, 
are shown in Figure 10. In general, equilibrium activities of alcohols and carboxylic acids ate 
quite high, especially at low temperatures. At temperatures relevant to hydrothermal systems, 
production of these simple organic compounds from alkyl thiols is quint significant. By comparing 
figures 10A and IOC with 10B and 10D, respectively, we find that the higher equilibrium activities 
of thiols formed from CO rather CO, have resulted in higher activities of alcohols and carboxylic 
acids. By comparing figure 10A to 10C and 10B to 10D, we can begin to understand how 
geochemical constraints can control organic compound production in hydrothermal systems. 
Higher activities of methanol, ethanol and acetic acid can be expected in Archean hydrothermal 
systems than in modem examples because of the different amounts of thiol, H, and H,S available 
in each system. Furthermore, in modem hydrothermal systems, methanol takes over as the most 
abundant organic compound of the three at very high temperatures, while in early Earth systems 
acetic acid is most abundant over the entire temperature range. Of course, these calculanons 
provide upper limits on the production of these compounds, as they do not take into account the 
formation of other organic compounds and the competition among all organic compounds for the 

carbon available in the system. 

Experimental Analogs 

While thiols have not been measured previously from hydrothermal systems, some 
experimental work may indicate that they could be intermediate species in such environments. 
Seewald (1994) demonstrated experimentally that predicted metastable equilbrium states do exist 
among organic species (aqueous hydrocarbons) under hydrothermal conditions. Because Seewald 
(1994) used sulfide-bearing minerals to control the oxidation state in his experiments, we believe 
that thiols were likely to have formed during his experiments and here we calculate the 
concentranons that should have been achieved based on his results. At 350 bars pressure and 
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■ f H S lanf was 6 1 ■ x 10’ molal and the H,(aq) concentration 
q25°C the reported concentrahon of H 2 S (aq) 

c n ff.™„ hv the nyrite-pyrrhodte-magnetite mineral 

was 0.36 x 10’ molal, consistent with buffeting by th py P 
assemblage used in the experiments. The ethane concentration was reponed to 
From the logarithm of the equilibrium constant for the reaction 


C 2 H 6 (aq) + H 2 S(aq) = C 2 H 5 SH(aq) + H 2 (aq), 

ethanethiol 


(ID 


6 

ethane 


r rh np H S and H we calculate that ethanethiol should 
Jong with the reported concentrations of ethane, H 2 S and 2 , 

lon * F , o Such a concentration 

r i n v 1 n-6 molal or approximately 2 ppm. 
lave reached a concentration of 1 9 x 10 molal, or PP 

ihould be detectable and would provide a direct test of these estimated data, ^ ^ 

achievement of metiistiible equilibrium among organic sulfur aud other organic ^ ^ ^ 

any other hydrothermal experiments conducted using sulfide minerals as re ox 
potential to genera, thiols that may be ptesen, in me, as, able equilibrium concentiations. 

CONCLUDING REMARKS 

Using previously developed methods, as well as some new collations 
have estimated standard partial molal thermodynamic properties of aqueous “ ' 

— - .so , 0 W calculation of reaction properties iuvolvmg — 
temperatures and pressures using be revised HKP equations of state. We n P 
incorporating these compounds into ge~hemical models of hydroma, organic geochenuc, 

"e results shown here imply dia, aqueous thiols may play a hey role in the orgamc 
geochemistry of hydrothermal systems. In general, the formauon of longer cham ^ 
or CO , HtS and H, is mom favorable ban . formation - — in 

h,w 250=C. This result lends support to previous model- 


hydrothermal systems that conclude drat the longer chain acids are more favorable to form at lower 
temperatures (Shock and Schulte, 1998). In addition, the formation of these compounds from 
carbon monoxide is mote facile that formation from carbon dioxide. For example, the acttvtty of 
methanethiol can theoretically reach values as high as 0.1 and 10 ’ from CO and CO„ respectively, 
in modem hydrothermal environments (see Figure 7). Extrapolating to the early Earth, the 
calculated values of methanethiol activities rise to as much as 1000 (from CO) and 10’ ! (from CO,). 
These high activities of methanethiol would result in a concomitant increase in the potential for 
production of organic compounds for which thiols serve as reaction intermediates. 

The formation of alcohols and carboxylic acids from Ms is thermodynamically favored, 
indicating that thiols may serve as reaction intermediates between CO or CO, and other organic 
compounds. The formation of organic compounds from thiols would create a positive feedback, tn 
which the consumption of thiols during organic compound synthesis would drop the thiol acttvtty 
below equilibrium values, thus driving further thiol synthesis and allowing thiols to persist in 
hydrothermal environments. We would encourage researchers making measurements of fluid 
samples from these systems to mclude these and other kinds of organic sulfur compounds in their 

analyses. 

While the formation of thiols from either CO or CO, is favored at temperatures and 
pressures relevant in hydrothermal systems, the amount of thiols that may be produced in these 
environments depends on the geochemical constraints imposed by each system. Thus, not only is 
the carbon source a major factor in determining tire equilibrium thiol activity, but the composition 
of die inorganic components of the system plays a role in tire distribution of organic compounds as 
well. For example, over a range of temperatures relevant to hydrothermal systems, the production 
of methanethiol from CO in an early Earth system is much greater than that in a modem system 
(see Figures 7 and 8). The importance of these geochemical constraints is evident when exploring 
the fotmation of alcohols and carboxylic acids using thiols as reaction intermediates. For example, 
in Archean hydrothermal systems, acetic acid is always more abundant titan ethanol or methanol, 
whereas in modem systems, methanol is the most abundant of these three above 300 C. 
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modem systems, the acdvity ratio of medranethiol to «dc add is . a, *PP— ly 250"C and 
the reaction to form acetic acid is thermodynamically favored. In early Earth systems, t e tger 
values of thiol activity would result in a concomitant increase m the potennal for production 

organic compounds for which thiols serve as reaction intermediates. 

The elevated potential of thiol syndesis in ear., Eardt hydrothermal systems places thtols ,n 
a hey role of organic compound synthesis. The presence of sulhrr in many essential biomolecules. 
including coenzyme-M, which is itself a thiol, supports this notion. Organic compound syn 
from thiols may have led to the first •ge.hem.cal metabolie cycles' from which mortem 
metabolism is derived. If life started in a hydrothermal environment, organic react, ons may avp 
taken place on sulfur-tearing mineral surfaces, which acted not only as catalysts, but a so as 
source for the sulfur that drove the first primitive metabolic cycles. 
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appendix 


Standard State Conventions 

The standard state employed to this study lot gases is urn. fugacity of the hypothetical W 
gas a, any temperature and 1 bar. The standard state convention for aqueous spectes adopte nr 
presen, study is one of uni. activity in a hypothetica. one moial solution referenced to tn ,n, e 
dilution at any pressure and temperature. The standard suite for H 2 0 calls for unit activity of th 

pure solvent at any pressure and temperature. 


Summary of the Revised HKF Equations of State 


One approach to calculating standard state properties for aqueous organ, c acids at ^ E 
pressures and temperatures which allows estimate and predictions for many spec.es invo ves e 
reused Helgeson-Kirkham-Flowers (HKF) equations of state (Helgeson et al., 1981; Tanger and 
Helgeson, 1988; Shock et ah. 1992). Data and pmameters for several groups of organic 
compounds are given by Shock and Helgeson <1990,. Schulte <1992,. 

21X11 i and Dale et ah (1997). The data and parameters for monocarboxyl.c an 

“Landes, e„ded.a„d combined wi, dam and parame te rs for d, carhop acids, hy,^ 

organic acids were used in tire thermodynamic calculations described ,n die text. 

The standard partial moial properties of aqueous species are expresse in ere 

equations of state in terms of structural and solvation contributions. A general statement o . 
summation for any standard partial moM property of the;,* aqueous species, Ej , can be expresse 


where the subscripts n and s refer to the nonsolvation (or structural) and solvation contributions, 
respectively. The revised HKF equations for the standard partial molal volume, V , and heat 
capacity, C;, are given by (Tanger and Helgeson, 1988) 

V° = AV n + AV° 
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where 5. «, a„ a,, a„ a„c„ c, and co represent species-dependent equation of - • * 

and 9 designate soivent-dependent parameters equ* to 2600 bars and 228 K for H A T and P 
stand for temperature and pressure, P, signifies the reference pressure of 1 bar, and Q, Y, and X 
indicate the Bom functions defined by 
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where £ stands for the dielectric constant of H 2 0. Equations (A-2) and (A- 3) have been used to 
regress experimental data for aqueous inorganic and organic electrolytes at elevated temperatures 
and pressures (Tanger and Heigeson, 1988; Shcck and Helgeson, 1988; 1990; Shock, . 

and pressures 1007- Prapaipong et al., 

Shock and Koretsky, 1993, 1995; Shock et al„ 1997; Sveqensky et al„ .997, 

1999 ), and are consistent with the approach of V and C p of aqueous electrolytes to 

critical point of H 2 0 (373.917°C and 220.46 bars). 



In the case of neutral aqueous species, the solvation terms in the revised HKF equations of 
state can be expressed in terms of the effecdve Bom coefficient (to,), which is given for the nth 

neutral species by (Shock et al., 1989) 

_ (A-7) 


where r„ and Z„. refer to the effective electrostattc radius and effecdve charge of the n,h neutral 
aqueousspecies, respectively, and „ = 1.66027 x 10 s (A cal mo. '). Regression of experimental 


data is consistent with the assumption that the effecdve Bom coefficients of neutral aqueous 
species are independent of temperature and pressure (Shock et al„ 1989; 1992, Shock and 
Helgeson, 1990; Svetjensky et al.. 1997; Shock, 1995; Schulte and Shock, 2001; Schulte etal., 
2001). I. follows that the revised HKF expressions for the solvation condibuiions to the standard 
partial molal volumes and heat capacities of neutral aqueous species are given by 
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As a consequence, the revised HKF expressions for die standard partial molal volume and heat 
capacity of neutral aqueous species correspond to (Shock et al., 1989) 
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where T indicates the reference temperature of 298. 1 5 K. 

^ apparent standi part, a, mo, a, Gibbs fee energy and enthalpy of formation are P ven 

by (Benson, 1976; Helgeson etal.. 1981; Tanger and Helgeson, 1988) 

(A-13) 
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where AG} and AH} stand for the 
enthalpy of formation at the reference conditions 
parenthetical terms represent the difference between 


of 1 bar and 298.15 K, respectively, and the 
the standard partial molal Gibbs free energy 



enthalpy of formation at the reference condittons and those a, the temperature and pressure of 
interest. I. follows from Eqns. (A-2), (A-3), (A-13), and (A-14) together with 
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and the partial isobaric temperature derivative of Eqn. (A-2), that we can write 
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Equations (A- 10), (A-ll). (A-17) and (A-18, are used to regress experiment^ d^^or^^eous 

neutral species at elevated temperatures and pressures, as well as to make pre ictio 

w rchnck et al. 1989; Shock and Helgeson, 

species for which experimental data are not aval a ’ , 2001). 

1990; Schulte and Shock, 1993; 2001; Shock, 1995; Dale et al., 1997; c 6 ” 

These calculations are consistent with the approach of C°, S\ V°, and A H/ of aqueous neutr 


species to « at the critical point of H 2 0. 


Table 1 . Summary of standard partial molal thermodynamic properties 
of aqueous alkyl thiols at 25°C and 1 bar, calculated with data taken 

from the literature. 



AG]* 

ah;* 

S° b 

methanethiol 

-1840.° 

-11650/ 

38.60 e 

ethanethiol 

-1100. c 

-17820/ 

45.50 c 

propanethiol 

1380/ 

-23320/ 

53.90 e 

butanethiol 

3560. c 

-28630/ 

61.30 e 

pentanethiol 

6150/ 




a. cal mol-' b. cal mol' 1 K l c. Calculated from AG° f (g) from Alberty et 
al. (1987) and AG; (goaq) from Abraham etal. (1990) d. Calculated 
from AH; (g) from Stull et al. (1969) and AHj (g^aq) from Abraham 
etal (1990) e. Calculated from AH; and AG] in the table, together 
with S° of the elements from Cox et al. (1989) f. Calculated from 
AG}(g) from Alberty et al. (1987) and AG r ° (g^aq) from Amoore and 

Buttery (1978) 
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Table 3. Reactions describing mineral assemblages that can buffer 
the activities of chemical species in hydrothermal environments. 



(Pyrrhotite-Pyrite-Magnetite) 


1.5 FeS + H 2 0 = 0.75 FeS 2 + 0.25 Fe 3 0 4 + H 2 
pyrrhotite P>' rite magnetite 

0.25 FeS 2 + 0.5 FeS + H 2 0 = 0.25 Fe 3 0 4 + H 2 S 
pyrite pyrrhotite magnetite 


FMQ 

(Fayalite-Magnetite-Quartz) 

1.5 Fe 2 Si0 4 + H 2 0 = Fe 3 0 4 + 1.5 Si0 2 + H 2 
fayalite magnetite quartz 


PM 

(Pyrrhotite-Magnetite) 


3 FeS + 4 H 2 0 = Fe 3 0 4 + 3 H 2 S + H 2 
pyrrhotite magnetite 


Table 4. Reactions forming methanol from methanethiol and ethanol and 
from ethanethio,, along with the reduced laws of mass action required 
logarithms of the equilibrium activities of methanol, ethano, and acettc acrd. 
resulting values of the log equilibrium activittes are plotted in Figure 



Mpfhanol F ormation 


CH 3 SH + H 2 0 = CH 3 OH + h 2 s 

methanethiol methanol 



Ffhannl Fo rmation 


c 2 h 5 sh + h 2 o = c 2 h 3 oh + h 2 s 

ethanethiol ethanol 


log a C2 H 5 OH 


log K + log &c 2 H $ SH ' a » 2 s 



Acid F ormation 

C 0 H 5 SH + 2H 2 0 = CH 3 COOH + H 2 S + 2H 2 
ethanethiol acetic acid 


lOg a CH 3 COOH 


log K - log a HjS - 2 log a Hj 


+ log a C2 H jS H 



Figure Captions 


Figure 1 Standard partial molal Gibbs free energy of formation of aqueous alkyl thiols as a 
function of the number of carbons in the alkyl chain. The free energy values for the aqueous 
species were calculated from the standard Gibbs free energies of formation of the gas (Alberty et 
al., 1987) and the free energy of solubility reactions (Abraham et al., 1990); see Table 2 for 

details. 

Figure 2. Standard partial molal enthalpy of formation of aqueous alkyl thiols as a function of the 
number of carbons in the alkyl chain. The enthalpy values for the aqueous species were calculated 
from the standard enthalpies of formation of the gas (Alberty et ah, 1987) and the free energy o 
solubility reactions (Abraham et ah, 1990); see Table 2 for details. 

Figure 3. Standard partial molal entropy of aqueous alkyl thiols as a function of the number of 
carbons in the alkyl chain. See Table 2 for details of how these values were calculated. 

Figure 4. Estimated values of the standard partial molal heat capacity (left-hand y-axis; open 
circles) and volume (right-hand y-axis; open squares) of aqueous alkyl thiols as a function of the 
number of carbons in the alkyl chain. Also plotted are the experimentally determined values o 
these properties for H 2 S (Shock et ah, 1992), which plot at n = 0 (filled symbols). 

Figure 5. Logarithm of equilibrium constants for reactions to form aqueous alkyl thiols from CO, 

H 2 S and H 2 as a function of temperature (at saturation pressure). The general reaction is shown 

for reference. 

Figure 6. Logarithm of equilibrium constants for reactions to form aqueous alkyl thiols from C0 2 , 
H 2 S and H 2 as a function of temperature (at saturation pressure). The general reaction is shown 
for reference. The letter “n,” corresponding to the number on the curves, refers to the number o 
carbons in the alkyl chain and determines the stoichiometric reaction coefficients for various 
species in the reactions (e.g. T refers to methanethiol, '2' to ethanethiol, etc.). 

Figure 7 Logarithm of methanethiol activity as a function of temperature for modem hydrothermal 
vent systems. The activity of H 2 S and the hydrogen fugacity are set by values in equilibrium with 
the PPM mineral assemblage. The concentration of C0 2 (13.8 mmol) for modem systems was 
taken from the literature for 1 1-13°N EPR vent fluids (Von Damm, 1990) and the CO concentration 

was estimated at 10’ 6 m. 


I 



Figure 8. Logarithm of me.hane.hiol vrith the FMQ 

Rsure 9 Logarithm of equilibrium coustants for reactions in Table 4, forming methanol from 
methanethiof and ethano, and acetic acid from ethanethtol at 250 bars pressure. 

and ethanol and acetic acid from ethanethiol as un assuming that the methanethiol and 

These values were calculated for modem hydro erm sy assuming that the 

, - nr . m Values for modem hydrothermal systems assuu & 
ethanethiol formed from CO. ) Q These values were calculated for early 

methanethiol and ethanethiol were formed from CO^ ^^^olfo^frornCO. D) 

Earth hydrothermal systems assuming the met ane 10 d ethanet hiol formed 

Values for early Earth hydrothermal systems assuming the methanethiol 

from C0 2 . 
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